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Abstract 

Yield  stresses,  allowable  stresses,  moment  capacities  (plastic  moments), 
external  loadings,  manufacturing  errors,...  are  not  given  fixed  quantities 
in  practice,  but  have  to  be  modelled  as  random  variables  with  a  certain 
joint  probability  distribution.  Hence,  problems  from  limit  (collapse)  load 
analysis  or  plastic  analysis  and  from  plastic  and  elastic  design  of  struc¬ 
tures  are  treated  in  the  framework  of  stochastic  optimization.  Using  espe¬ 
cially  reliability-oriented  optimization  methods,  the  behavioral  constraints 
are  quantified  by  means  of  the  corresponding  probability  p^  of  survival. 
Lower  bounds  for  p^  are  obtained  by  selecting  certain  redundants  in  the 
vector  of  internal  forces/bending-moments;  moreover,  upper  bounds  for  p^ 
are  constructed  by  considering  a  pair  of  dual  linear  proprams  for  the  op- 
timizational  representation  of  the  yield  or  safety  constraints .  Whereas  the 
probability  p^  can  be  computed  e.g.  by  sampling  methods  or  by  asymptotic 
expansion  techniques  based  on  Laplace  integral  representations  of  certain 
multiple  integrals,  efficient  techniques  for  the  computation  of  the  sen¬ 
sitivities  (of  various  orders)  of  p^  with  respect  to  input  or  design  vari¬ 
ables  X  and  load  factors  A  have  yet  to  be  developed.  Hence,  several  new 
techniques  (e.g.  Transformation  Method,  Stochastic  Completion  Technique) 
are  suggested  for  the  numerical  computation  of  derivatives  of  p^  with  re¬ 
spect  to  (X, A) . 
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1.  Limit  (collapse)  load  analysis  of  structures  as  a  linear  program¬ 
ming  problem. 

The  collapse  load  can  be  defined  [3],  [5], [18]  "as  the  load  required  to 
generate  enough  number  of  local  plastic  yield  points  (referred  as  plastic 
hinges  for  bending  type  members)  to  cause  the  structure  to  become  a  mech¬ 
anism  and  develop  excessive  deflections".  Assuming  that  the  material  be¬ 
haves  as  an  elastic-perfectly  plastic  material  [4], [15],  a  conservative 
estimate  of  the  collapse  load  factor  is  based  on  the  following  formu¬ 
lation  as  a  linear  program  (LP) : 

maximize  A  (1) 

s .  t . 


F^  <  F  <  F^ 
CF  =  AR  . 


(1.1) 

(1.2) 


Here,  (1.2,)  is  the  equilibrium  equation  of  a  statically  indeterminate 
loaded  st;ru6ture  involving  an  mxn  matrix  C  -=  (c^j),  m<n,  of  given  coeffi¬ 
cients  c_ ,  l<i<m,  l<j<n,  depending  on  the  undeformed  geometry  of  the 
structure  having  n^  members  (elements) ;  we  suppose  that  rankC^m.  Further¬ 
more,  R  is  an  external  load  m- vector,  and  F  denotes  the  n- vector  of  in- 
o 

ternal  forces  and  bending -moments  in  the  relevant  points  (sections,  nodes 
or  elements)  of  the  structure.  Finally,  (1.1)  are  the  yield  conditions 
with  the  vector  of  lower  and  upper  bounds  F^,F^. 

For  a  plane  or  spatial  truss  [6],  [16]  we  have  that  n=n^,  the  matrix  C 
contains  the  direction  cosines  of  the  members,  and  F  involves  only  the 
normal  (axial)  forces,  moreover, 

F^:=<t^.A.,  FV;=  a^.A.,  j-1 . n(=n  ),  (1.3) 

J  yj  J  J  yj  J  ^  o^’  ^  y 

where  A.  is  the  (given)  cross-sectional  area,  and  a  a  resp.,  denotes 
the  yield  stress  in  compression  (negative  values)  and  tension  (positive 
values)  of  the  j-th  member  of  the  truss.  In  case  of  a  plane  frame,  F  is 
composed  of  subvectors  [16] 


/k«\ 
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(1.4) 
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where  denotes  the  normal  (axial)  force,  and 

the  bending-moments  at  the  positive,  negative  end  of  the  k-th  member.  In 
this  case  F^,  contain  -  for  each  member  k  -  the  subvectors 


(1.4.1) 


resp.  ,  where  k=«l,...,n^,  denote  the  plastic  moments  (moment  capac¬ 

ities)  given  [4], [15]  by 

and  “  ^kpl^^^  plastic  section  modulus  of  the  cross-section  of 

the  k-th  member  (beam)  with  respect  to  the  local  z-axis. 

For  a  spatial  frame  [6], [16],  corresponding  to  the  k-th  member  (beam),  F 
contains  the  subvector 

where  t^  is  the  normal  (axial)  force,  m^^  the  twisting  moment,  and  m^^, 
m^- ,m^^,mj^-  denote  four  bending  moments  with  respect  to  the  local  y-,z- 
axis  at  the  positive,  negative  end  of  the  beam,  respectively.  Finally, 
corresponding  to  (1 . 3) , (1 .4) - (1 .4 . 2) ,  the  bounds  F^,F^  for  F  are  given  by 


(1.5) 


^  “  ^‘"yk\’^pl’^pl’\pl’^pl’\pl^' ’ 


(1.5.1) 


(1.5.2) 


where,  cf .  [4] , [15] , 

^pl‘”  ’’yk'^^pl’  ^pl'”  °'yk'*^^pl’  \pl'“  ‘"yk\pl  (1.5.3) 

are  the  plastic  moments  of  the  cross-section  of  the  k-th  element  with  re¬ 
spect  to  the  local  twisting  axis,  the  local  y-,z-axis,  respectively.  In 


(1.5.3),  wPp^=wPp^(X)  and  P'” 

lar,  axial  modulus  of  the  cross-sectional  area  of  the  k-th  beam  and  r 

U  ^ 

denotes  the  yield  stress  with  respect  to  torsion;  we  suppose  r  ,=cr  ,  . 


yk  '^yk* 


Remark  1.1 


Possible  plastic  hinges  [4], [7], [15]  are  taken  into  account  by  insert¬ 
ing  appropriate  eccentricities  »  k=l,...,n^,  with  » 

where  is  the  length  of  the  k-th  beam. 

Remark  1.2 

Working  with  more  general  yield  polygons  [1] , [17] , [18] ,  the  stress  con¬ 
dition  (1.1)  is  replaced  by  the  more  general  system  of  inequalities 

H(F|!)‘S<h.  (1.1a) 

^  u  u 

Here,  (H,h)  is  a  given  i/x(n+l)  matrix,  and  F^:  =  (Fj5^j)  denotes  the  nxn 
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diagonal  matrix  of  principal  axial  and  bending  plastic  capacities 

r 

VA, 


P":. 

J 


F^: 

J 


W,Jp. 


yk. 

J  J  J  J 

where  indices  as  arising  in  (1.4.1)-(1.5.3).  The  more  general 

case  (1.1a)  can  be  treated  by  similar  methods  as  the  case  (1.1)  which  is 
considered  here. 


2.  Plastic  and  elastic  design  of  structures 

In  the  plastic  design  of  trusses,  frames  [5]  having  n^  members ,  the  n- 
vectors  F^,F^  of  lower  and  upper  bounds 


-  F^(a^,a^,X),  F^  =  F^Ca^.a'^.X) 

y  y  y  y 


(2) 


for  the  n-vector  F  of  internal  member  forces  and  bending-moments  F.,  j=l,. 
. . ,n,  are  determined  [3]  ,  [5]  by  the  yield  stresses,  i.e.  compressive  lim¬ 
iting  stresses  (negative  values)  =  (a^-,...,a^  )',  the  tensile  yield 

y  y 

u  ,  u  u  ,  ,  ,  , 

stresses  a  ==(c7-,...,a  )',  and  the  r-vector 

y  yi  y^Q 

X  =  (X^,X2,...,X^)'  (2.1) 

of  design  variables  of  the  structure.  In  case  of  trusses  we  simply  have 
that ,  cf .  (1.3), 

F^  =  <^ydA(X)  “  A(xyj;,  F^  =  -  A(X)^aj;.  (2.2) 

where  n=n^,  and  cr^^,cr^^  denote  the  nxn  diagonal  matrices  having  the  diago¬ 
nal  elements  resp.  ,  i-“l,...,n:  moreover, 

yj  yj  *  J  » 


A(X)  «  (A^(X),...,A^(X))' 


(2.3) 


is  the  n-vector  of  cross-sectional  areas  Aj*Aj(X),  j"»l,.,,,n,  depending  on 


the 


r-vector  X  of  design  variables  X^,  k*=l,,..,r,  and  A(X)^  denotes  the 


nxn  diagonal  matrix  having  the  diagonal  elements  Aj*=Aj(X)  ,  l<j<n. 

Corresponding  to’ (1.2),  here  we  have  again  the  equilibrium  equation 
CF  =  R^,  (2.4) 

where  describes  [5]  the  ultimate  load  (representing  constant  external 
loads  or  self-weight  expressed  in  linear  terms  of  A(X)). 

The  plastic  design  of  structures  can  be  represented  then  [1],[2],[5]  by 
the  following  optimization  problem 

min  G(X)  (3) 

s .  t . 

T  T  TT  TT  T  TT 

(3.1) 


F^(aJ;.aj;.X)  <  F  <  FU(aJ;,aU.X) 


CF  =  R  , 
u 


(3.2) 

where  G=G(X)  is  a  certain  objective  function,  e.g.  the  volume  or  weight  of 
the  structure. 
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Remark  2.1 

As  mentioned  in  Remark  1.2,  working  with  more  general  yield  polygons, 
(3.1)  is  replaced  by  the  condition 


H(F^(ay,X)^)'^F  <  h. 


(3.1a) 


For  the  elastic  design  we  have  to  replace  the  yield  stresses  0^,0^  by 

L  u  .  .  y  y 

the  allowable  stresses  o  .a  y  and  instead  of  ultimate  loads  we  consider 

a  a 

service  loads  R^.  Hence,  instead  of  (3)  we  get  the  related  program 
min  G(X) 


F^(a^,(7^,X)  <  F  <  F^Cff^.a^.X) 

Si  a  Si  cL 


(4.1) 


X^  ^  X  S  X^ 

where  X^,X^  still  denote  lower  and  upper  bounds  for  X. 


(4.3) 


3.  Analysis  and  design  of  structures  in  case  of  random  data 
In  practice,  yield  stresses,  allowable  stresses,  the  loads  applied  to 
the  structure,  other  material  properties  and  the  manufacturing  errors  are 
not  given  fixed  quantities,  but  must  be  treated  as  random  variables  on  a 
certain  probability  space  (0,W,P).  Hence,  (1),(3),(4)  are  stochastic  lin¬ 
ear/  nonlinear  programs  (SLP/SNLP)  which  obviously  have  the  same  basic 


structure  represented  by  a  random  objective  function 

Z(w)  G(w,X)  ,  a;  G  n,  (5) 

and  by  stochastic  constraints  of  the  type 

CF  »  R(w)  (6) 

F^(w)  <  F  <  f''^(w)  ,  (7) 

where  R=R(w)  ,  o)  e  fi,  is  a  random  load  m-vector  given  by 
R(w)  =  AR^(w),  R(a;)  =  f  R(^)  “ 

resp.,  and  for  the  n-vector  F=(Fj)  of  internal  member  forces  and  bending- 
moments  we  have  the  n- vectors  of  random  bounds 

F^(w)  =  F^(a>,X),  F^(w)  -  F^(c<j,X),  w  e  0,  (7.1) 

depending  on  an  r-vector  X  of  design  variables  X,  ,  k=l,...,r. 


Obviously,  each  realization  of  the  random  element  cogO  yields  new  loading 
conditions,  represented  by  the  vector  R=R(a>) ,  and  therefore,  cf.  (6),  new 
arrangements  F=F(w)  of  internal  member  forces  and  bending-moments.  Hence, 
in  the  present  case  of  "multiple  loadings"  caused  by  the  random  variations 
of  R=R(c(;)  ,  the  survival  of  the  structure,  i.e.  the  existence  of  certain 
arrangements  F(o))  of  internal  member  forces/bending-moments  not  overwhelm¬ 
ing  the  strength  of  the  structure,  can  be  evaluated  by  the  probability  of 
survival 
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(8.2) 


(8.3) 


p  P(There  is  F“F(w)  such  that  CF(w)  -  R(w) 

^  L  U 

and  F  (cj)  <  F(o))  <  F  (w) )  , 

see  [2], [14],  assuming  that,  cf.  (26) , (33) , 

S(F^( . )  »F^(  • )  |R(  • ))  :=  {c.)en:  there  is  a  vector  F=F(w) 

fulfilling  (6)  and  (7)) 

is  a  measurable  set.  Denoting  by  [F,^,F^]  the  n-dimensional  interval 

[F^.F^]:-  {Fsk”;  F^<F<F^) ,  (8.2) 

we  find  that 

S(F^(.),F^(.),R(.))  -  {toGft:  R(a))  e  C[F^(w)  ,F^(w)  ] )  (8.3) 

with  C[F^,F^]  =  (CF:  F^<F<F^) ,  and  therefore 

-  P(S(F^(.),F^(.),R(.)))  =  P(R(a))GC[F\w),F^((^)])  (9) 

=  /P(R(a))GC[F^,F^]  |F^((j)=F^,F^(o)=F'^)7r(dF^,dF^)  , 
where  tt  denotes  the  distribution  of  the  bounds  (F^(w) , F^(w) ) .  Since  the 
bounds  F^,F^  in  (7)  depend  also  on  the  vector  X  of  design  variables  X^, 

k=l . r,  cf.  (7.1),  we  have  p^=»P(X)  with  the  probability  function 

P(X)  =  P(R(w)  e  C[F^(w,X),F^(w,X)]);  (10) 

furthermore,  if  the  external  load  R=R(w)  is  given  by 

R((^)  -  R(a>,A):  -  S  A  R^^^co)  ,  A:-  (A  . A  )'.  (10.1) 

1=1,7)  ^  "*R 

with  random  m-vectors  R^  R^  t  i=l ,  .  .  .  ,mj^,  and  deterministic  coeffi¬ 

cients  A^,  i=l ,  .  .  .  ,mj^,  then  p^=P(X,A)  with 
^  ™R  .  ^ 

P(A,X):  =  P(  S  A .R^^\w)  e  C[F^(w,X),f'^(w,X)]).  (10.2) 

i=l 

Especially,  in  case  of  trusses,  see  (1) , (3) , (4)  and  (2.2),  for  dealing 
with  the  probability  of  survival  p^  we  have  the  following  probability 
functions: 


(10.2) 


P^(A,X) 

Pu(X) 

Ps(X) 


-  P(AR^(w)  G  C[A(X)^ay(w).A(X)^/(w)]),  AgK 
=  P(R^(a))  G  C[A(X)^ay(u.),A(X)^<7y(a))]) 

=  P(Rg(w)  G  C[A(X)^a^(w),A(X)^ff^(«)]). 


(10.3) 

(10.4) 


(10.5) 


Since  p^*P(X)  are  very  complicated  expressions  in  general,  in  the  fol¬ 
lowing  we  are  looking  for  approximations  (lower,  upper  bounds)  of  p^“P(X) 
by  simpler  probability  functions. 


4.  Lower  and  upper  bounds  for  p^  =  P(X) 

According  to  (8.3)  we  have  that 

T  U  ^  T  II 

S(f'^(.),f‘'(.).R(.))C  r>S.(F^.).F''(.)). 

i-1  ^ 


where 


S^(F^(.),F^(.),R(.)):-  {wed:  R^(o})  e  C^[F^(w)  ,f'^(w)]  ) ,  i-1 . m,(ll.l) 

and  C.  denotes  the  i-th  row  of  C.  Hence,  we  get  the  inequality 


P(X)  <  min  P.(X), 
l<i<m  ^ 


where 


P^(X):-  P(Rj^(w)  e  C^(F^w,X),F^(u,.X)]) 


(12.1) 


Note  that  related  inequalities  (with  more  exact  bounds)  follow  from  more 
general  Bonferroni- type  inequalities  [10] .  We  find  that 

C^[FV,X),F^(a,,X)]  =  [7^(w,X),7V(w,X)]  (12.2) 

is  an  interval  in  having  the  bounds 


7^(w,X)  =»  rain  C.G*'(w,X) 

^  l<t^  ^ 

7V(w,X)  *=  max  C.G^(a),X), 
^  l<t<J  ^ 


(12.3) 


where  G^-'G^  (a>,X)  ,  i=l,...,J,  are  the  extreme  points  of  the  interval 

[F^(w,X) ,F^(w,X) ] .  Since  the  components  Gt(w,X),  j«l,2,...,n,  of  G^(w,X), 

L  J  U 

t“l,...,J,  are  certain  elements  of  (F  (co,X),F  (w,X)),  the  measurability  of 
the  bounds  F^,F^  with  respect  to  wG(n,0t,P)  yields  the  measurability  of  7^ 
and  7^,  i=-l,...,m,  with  respect  to  w.  Hence,  (F^( . ) , F^( . ) , R( . ) ) ,  i=l,.. 

. . ,m,  are  measurable  sets,  and  we  may  write 

P^(X)  -  P(7^(w,X)  <  Rj^(w)  <  7^(0), X)),  i-l,...,in.  (12.4) 


4.1.  Lower  bounds  by  selection  of  redundants 

For  the  construction  of  lower  bounds  for  P(X)  ,  the  vector  F«F(ci))  is 
partitioned 


into  a  certain  (n-m) -vector  N=(F.  )  of  redundants  F  ,  i=l,...,n-m, 

l<jg<n-m  ^Jl 

and  an  m-vector  F^  of  statically  determined  member  forces/bending-moments. 
Hence,  with  a  corresponding  partition  of  the  mxn  matrix  C  into  mxm,  mx(n-m) 


submatrices  Cj,C^^,  resp.,  where 
rankC^  =  rankC  =  m, 

the  equilibrium  equation  (6)  yields  for  F(co)  the  representation 


(13.1) 


C*^R(w) 


Consequently,  selecting  for  each  ci>gQ  a  vector  of  redundants 

N=N(a))  =  (F.  (cj))  such  that  N(.)  is  a  measurable  function  on  P), 

l<i<n-m 

we  get 

S(F^(..X),F^(.,X),R(.))DS(X,N(.),R(.)),  (IS] 

where  S(X,N( . ) ,R( . ) )  is  the  measurable  set  given  by 


C  "(R((o)-C  N(w)) 

S(X,N(.),R(.)):=  {coed:  F^(w.X)  <  ( 


)  <  F^(w.X)).  (15.1) 


Thus,  for  P(X)  we  find  the  lower  bound 
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P(X)  >  P(X,N(.)), 
where 


P(X,N(.)):=  P 


Fj(a,,X)  <  C'^(R(w)  -  CjjN(w))  <  Fj(w.X) 
Fjj(a>,X  )  <  N(a>)  <  Fjj(w,X) 


(16) 


(16.1) 


and  Fj.Fjj  and  Fj.Fjj  denotes  the  partition  of  F^.F^,  resp.  ,  correspond- 
ing  to  (13).  Note  that  the  inequality  (16)  holds  for  any  choice 

N=(F.  )  of  an  (n-m) -subvector  of  redundants  such  that  (13.1)  holds 

l<i<n-m 

and  any  representation  of  N  as  a  random  vector  N=N(a))  on  (n,(K,P);  espe¬ 
cially,  N  can  be  selected  as  a  deterministic  vector  of  redundants: 

N(w)  “  z  a.s.  (almost  sure),  (17) 

where  ze!R^  is  a  deterministic  vector;  in  this  case  we  set  P(X,N(.))  = 
P(X,z). 

4.1.1.  Special  cases,  a)  In  case  of  trusses,  cf.  (1.3), (2. 2),  we  have 
that 


Fj(w,X)  -  Aj(X) ^<7^(10),  Fj(w.X)  =  Aj(X)^aj(u,) 

Fjj(w,X)  -  Ajj.(X)^c7jj.(a,),  Fjj((0.X)  -  Ajj(w)^ajj(co). 


(18) 


where  A^ .A^j .cr^j are  the  partitions  of  A.a^.a^,  resp.,  corre¬ 
sponding  to  the  partition  F^.F^^  of  F,  and  Aj.(X)^  denotes  the  diagonal 
matrix  which  has  the  components  of  Aj(X)  as  its  diagonal  elements.  Thus, 
we  get 


P(X,z) 


=  P 


Aj(X)^c7j(w)  <  C'^(R(w)  -  CjjZ)  <  Aj(X)^aj(cc) 
Ajj(X)^(7jj(w)  <  z  i  Aj.j(X)^ffj(w) 


a^(a))  <  Aj(X)‘^C'^(R(w)  -  CjjZ)  <  ^^(w) 


(18.1) 


b)  Suppose  that  the  partition  of  F^,F^  into 

chosen  such  that  (F^(w,X) , F^ (w,X) ) ,  (F^ ^ (w) , F^ ^ (w) )  are  stochastically 
independent.  If  (17)  holds,  then 

P(X,z):=  P(Fj(w,X)  <  c‘^(R(w)  -  CjjZ)  S  Fj(w,X))  (18.2) 

X  P(Fjj(w,X)  <  z  <  Fjj(w,X)). 


5.  Failure  modes,  limit  state  functions  and  upper  bounds  for  P(X) 
According  to  (8), (10)  we  have  that 
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P(X)  =  P(There  is  F“F(a))  such  that  CF(a>)«R(a))  and 


Fj(<^)  -  Fj(w.X)  ^  0,  j-1, 


,n 


(19) 


F^(u,X)  -  F.{o>)  <  0, 

where  we  suppose  that  all  bounds  fV»fV  are  finite,  i.e. 

L  U  J  J 

-«><  F.(aj,X)  <  F.(w,X)  <-H=o  a.s.,  l<j<n,  for  all  X  under  consideration. 


(19.1) 


Defining 

t(w,F(a)),X):=  max  {F  (w) -F^(w,X)  ,  F^^.X) -F  (w) ) ,  (19.2) 

1<  j  <n  J  J  J  J 

we  obtain 

P(X)  -*  P(There  is  F=F(w)  such  that  CF(cj)=R(a>)  and  t (w, F(w)  ,X)<0) 

-  P(inf{t(w,F(w)  ,X)  :  CF(w)«R(a))  )  <  0)  (20) 

»  P(t'^(c^,X)  <  0), 


where 

t*(w,X):-  inf{t(a;,F(a>)  ,X)  :  CF(a))=R(tj)  ) 
is  the  minimal  value  of  the  program 
min  t(a>, F(u?)  ,X) 
s.  t. 

CF(w)  »  R(w) 

being  equivalent  to  the  linear  program 
minimize  t 


(20.1) 

(20.2) 


(21) 


F^-Fj(a),X)  -  t<0,  j-l,....n 
Fj(uj,X)  -  Fj  -  t<0,  j-l,...,n 
CF(w)  -  R(w) 


(21.1) 

(21.2) 

(21.3) 


with  the  variables  F-,F^,...,F  ,t. 

1  2  n 

Because  of  condition  (19.1),  for  each  (c»>,X)  we  get 

t(w,F(w),X)  >  max  i(FV(a),X) -fV(w,X)  )  >-»  (21.4) 

l<j<n  J  J 

for  arbitrary  F(co)  ;  hence,  the  objective  function  of  the  linear  program 
(21)  is  bounded  from  below  for  all  (cj,X).  Since  the  LP  (21)  always  has  a 

F 

feasible  solution,  for  each  (co,X)  an  optimal  solution  (  ^)  of  (21)  is 

t 


guaranteed,  and  we  have  that 

t*  -  t(w,F*(u>)  ,X)  -  t*(w,X).  (22) 

* 

Consequently,  by  means  of  duality  theory  the  optimal  value  t  (a),X)  of 
the  equivalent  programs  (20.2)  and  (21)  can  be  represented  also  by  the 
optimal  value  of  the  dual  program  of  (21)  given  by 
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max  R(w)'u  -  F^(w,X)'u'*’  +  F^(u),X)'u’ 


C'u  -  u  +  V  =  0 


(23.1) 

I'u"^  +  I'u'  -  1  (23.2) 

u^>0,v>0,  (23.3) 

where  is  not  restricted. 

Remark  5.1 

Obviously,  (23.1)  is  the  member -node (or  joint)  displacement  equation. 
According  to  its  mechanical  meaning,  we  call  (21) , (23) ,  resp. ,  the  static, 
kinematic  linear  program  (LP) ,  cf.  [2], [14]. 

Having  (23),  t  (ci),X)  reads 

/R(a))  Y  /u\ 

t*(a;,X)  =  max{  I -F^(w,X)  j  6:  I  u^j  ^:8  G  A  ),  (24) 

\F^a;,X)/  \v’/  ^ 

where  denotes  the  convex  polyhedron  in  K  represented  by  the  con¬ 
straints  (23.1)  -(23.  3)  of  the  LP  (23).  Taking  any  subset  A- C  A  of  A  ,  and 


I'-  o 


defining  then  t^(w,X)  by 

/  Y 

t*(«,X):=  sup{( -F^(w,X)  U;  «  e  A^), 
\f^(w.X)/ 

we  get 

t*(w,X)  >  t*(w,X), 

which  yields  for  P(X)  the  following  upper  bound: 
P(X)  ■=  P(t*(u,X)<0)  <  P(t*(w.X)  <  0). 
Moreover,  if 

.i=l . V 

\u-<'^>/ 


(25.1) 


(25.2) 


denote  the  extreme  points  of  the  convex  polyhedron  A^,  then 

t*(w,X)  =  max  R(w)'u^'*^  -  F^(w,X) 'u'*'^'*^  +  f’^(w,X) 'u' 
!<£<£ 


■st  L  U 

which  shows  that  t  ( .  ,X)  is  measurable.  Hence,  S(F  (.),F  (.),R(.))  =»  {cjeQ: 

t  (w,X)<0)  is  measurable,  cf.  (10) , (20) ,  and  we  get 

P(X)  -  P(R(w)'u^'*^  -  F^(w,X) +  F^(w,X) 'u"^'*^  <  0,  l<i<je^) .  (26.1) 
According  to  (6), (7),  the  survival,  failure,  resp.,  of  the  underlying 
structure  can  be  described  by  the  inequality 
k  k 

t  (w,X)  <  0,  t  (w,X)  >  0,  respectively. 

Thus,  the  structure  fails  if  and  only  if 


(27) 
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R(«)'u^'®^  +  -  F^(w,X)'G'^'*^  >  0  for  at  (27.1) 

least  one  l<i<>8  ; 

o 

obviously,  (27.1)  represents  the  different  failure  modes  of  the  structure. 

(i,) 

Having  a  certain  nviinber  of  basic  solutions  S  ,  r=l,...,i^,  of 

the  LP  (23),  and  defining 

*  U  L 

t  (co,X):=  max  R(a))'u  -  F^(o>,X)^n  +  f"^(w,X)'u  ,  (28) 

^  l<r<i^ 

corresponding  to  (25.1),  here  we  get 

t*(to,X)  >  t*(w,X)  (28.1) 

and  therefore 

P(X)  -  P(t*(w,X)  <  0)  S  P(t*(w,X)  <  0).  (28.2) 

6.  The  probability  of  failure  p^ 

According  to  (8) , (10) , (20)  and  (27),  for  the  probability  of  failure 
p^:=l-p^*=l-P(X)  we  obtain 
p^  =  P(t’^(a?,X)  >  0) 

-  P(R(w)'u  -  F^(u,X)'u'^  +  F^(cj,X)'u'  >  0  for  at 

=  P(R(w)'u^'*^  -  F;'^(w,X)'u'^^'^^  +  F^(«,X)'u’^'^^  >  0  for  at  least 

one  l<j2<i  ) 
o 

«  P(  U  F  (X)), 

>e=i 

where  F^(X)  denotes  the  ;?-th  failure  domain 

F^(X):=  (weO:  R(w)'u^'*^  -  F^(w,X) +  F^(a),X) 'u‘ >  0) 

-  (wen:  M^(w,X)  <  0) 

with  the  corresponding  limit  state  function 

M^(w,X):=  F^(w,X)'u'^^'*^  -  F^(w,X)'u'^'^^  -  R(w)'u^'*\ 

i=l,...,i^;  especially,  for  trusses,  cf.  (2.2),  (18.1),  we  find 

M^(w,X):-  (7^(w)'A(X)^u‘*'^'*^-  ct’^(w)'A(X)^u'^'*^  -  R(w)'u^'^\ 

Using  known  inequalities  for  probabilities,  for  p-  we  find  the  bounds 
i 

o 

max  p-  .  <  p.  <  S  p.  (30) 

l<f<i  ’  i=l  ' 

o 

where  p^  ^  is  given  by 

Pf,i*^  P(F^(X))  ==  P(M^(a>,X)  <  0) 

■=  P(F^(w.X) -  F^(w,X)'u'^'®^  <  R(w)'u^'^^ 

=  1  -  P(R(w)'u^'^^  i  f'^(w,X)'u‘'’^'*^  -  F^(w,X) 'u*^"^^)  , 


(29.1) 

(29.2) 

(29.3) 


w  \ 

least  one  I  u 


u 


roQ\ 


(30.1) 


-  11  - 

and  sharper  bounds  can  be  obtained  by  using  more  genral  Bonferroni- type 
inequalities  for  probabilities. 


7.  Representation  of  by  using  cones 

According  to  (9)  we  have  that 

=  P(R(a,)  G  C[F^(a.).F^(w)]), 

^  L  U 

where  [F  ,F  ]  is  given  by  (8.2).  Representing  therefore  the  vector  F  of 
internal  member  forces/bending-moments  by 
F  =  F^+AF^  =  f’^-AF^ 

with  n- vectors  AF^,AF^>0,  the  condition  R  e  C[F^,F^]  can  be  represented  by 


R  -  CF 


=  -  CAF 


-  (F^-Fb  =  -  AF^  -  AF^ 


AF^O,  AF^O. 


Thus,  we  consider  the  cone  Y  defined  by 


VI  I / \ AF 
2n 


2  y^:  k«l , . . . , 2n) , 

(k) 

where  the  cone -generators  y^  ,  k=l, . . . ,2n,  are  given  by 


,  l<k<n,  y^ 


1,  n<k<2n, 


(32.1) 


and  c.  ,  e.  denotes  the  k-th  column  of  C,  of  the  nxn  identity  matrix  I, 

L  U 

respectively.  Having  Y^,  the  set  S(F  (.),F  (.),R(.))  defined  in  (8.1)  can 
be  described  by 


S(F^(.),F^(.),R(,))  =  {wGa:(. 


R(a))-CF  (w,X)  \ 

IT  T  h 

-f''(£.),X)+fV,X)  J 


which  shows  again  the  measurability  of  S(F^( . ) , F^( . ) ,R( . ) ) . 

Moreover,  the  probability  function  P=P(A,X)  representing  (10)- (10.3)  can 
be  given  by 


P(A,X)  *  P( 


Remark  7.1 


CF  (w,X)  -  AR^(a>) 
F^Cw.X)  -  F^w.X) 


^  V- 


(33.1) 


The  cone  contains  all  pairs  (CAF^, AF^+AF^)  of  admissible  pairs  of 
external  load/strength  increments. 

According  to  the  representation  (32)  of  Y  ,  there  are  a  finite  number 

m+n  °  (£) 

of  boundary  hyperplanes  in  R™  represented  by  vectors  • 

i=l, . . . ,  such  that 

Yg  =  {y  =  (^)  G  w'w^'*^+  v'v^'^^>0,  l<i<.e^).  (34 

Hence,  because  of  (33)  and  (34),  the  survival  of  the  structure  can  be 


represented  also  by  the  inequalities 

(R(w)-CF^(w,X))'w^'^^+  (-F^(w,X)  +  F^(w,X))'v^'*^<  0, 
which  yields 

(R(w)'w^'*^-  F^(w,X)' (C'w^‘*^+v^'*^)  +  Fhw,X)'v^'*^<  0,  l<^<r,  ( 

and  therefore 

p^-  P(R(«)'w^'^^-  F^(w,X)'(C'w^'*^  +  +  F^(w,X)'v^'*^sO,  l<i<i^)  . 

(3; 

Obviously,  the  conditions  for  structural  safety  given  by  (27)  and  (35) 
coincide . 

(£) 

For  arbitrary  subset  Y  ,  i-1,2,  such  that 

O  o  O 
we  get 


where  the  bounds  ,  i-1,2,  are  defined  by 

■  ?<">  Uy<«V  1-1.2. 

®  \\F^(a.,X)  -  F^(w,X)  j  °  j 

7-1.  Construction  of  approxijnating  cones 

Suppose  next  to  that  we  have  a  cone  axis  (center  or  middle  line) 
g  {Ay:  A>0) 

generated  by  a  vector  y+^»  which  is  defined  later  on. 


^.(6) 


;{6) 


I 


^'7 
/  \A 


I 

L  - 


,(i)  \  /  // 

■  \  I'  /  /// 

\  I  /  vX 

Wlr 


/  /Y' 


Fig.  7.1.  Cone  Y  with  cone -generators  y'^  ,  i==l . i=l,...,6,  cone 

o  — 

axis  g,  hyperplane  E  and  convex  polyhedron  K  =conv(y  , . . . ,y 
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Let  denote  the  hyperplane 

(y-y)'y  -  o  <—>  y'y  -  I lyl  I' 


(37.1) 


through  y  and  orthogonal  to  axis  g.  Consider  then  the  points  y'"  'e  Y  , 
i-l,...,2n,  lying  on  E^  and  on  the  straight  lines  through  0  and  y  , 
i»l, , . . ,2n,  hence, 

-(«..ijill!  w  J.1  2„ 

y  / 

see  Figure  7.1.  Obviously,  we  get 

y^^^'y>0,  . 2n.  (3) 


(38.1) 


as  a  condition  for  y.  Since  the  equation  S  oc.y'^  -  0,  ot,>0,  k-l,...,2n, 

k-1  ^  ^ 

cf.  (32), (32.1),  has  no  nonzero  solution  a,  according  to  Gordan's  trans¬ 
position  theorem,  system  (38.1)  has  solutions  yfO.  Let 

^  ~(1)  ~(2)  ~(2n), 

K  conv{y"  ,y^  , . . . ,y^  ') 


denote  the  convex  polyhedron  oh  E^  generated  by  the  points  y'  '',i-l,.. 
.,2n.  By  (38.1)  and  (38.2)  we  get 


y'y>0  for  all  yeK^.  (38.1a) 

According  to  (38)  and  (38.1),  the  cone  can  be  represented  by 

Y  -  K  -  {Ay:  A>0,  yEK  ) .  (38.3) 

A>o 

Based  on  the  above  representation  of  Y^,  approximations  Y^  of  Y^  of  the 
type  (36)  can  be  obtained  by  replacing  the  generating  polyhedron  in  Y^ 
by  suitable  approximations,  e.g. ,  ellipsoides,  spheres,  denoted  by  K^: 


Fig.  7.2.  Convex  polyhedron  in  hyperplane  E^  and  (second  order)  ap- 
proximnations  \  j-1,2 
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a)  Ellipsoidal  approximations:  According  to  (38.1a),  an  outer  ellip¬ 
soidal  approximation  can  be  defined,  for  j-2,  by 

{yeE^.y.y>0,  (y-y)»r^j^r^J^'(y-y)  :S  1),  (39) 

(2) 

where  the  (n+m)x(n+m)  matrix  -r  is  chosen  such  that 
(y(^).y)'rr'(y^^^-y)  <  1,  i-l,...,2n 

y^y  -  0  — >  (ry)'y  -  0  (39.1) 

2  (1- (y^^^ -y) 'rr' (y^"^^ -y) )  - min. 

i-1 

Conversely,  an  "inner"  ellipsoidal  approximation  of  can  be  de¬ 
termined,  cf.  (39),  by  choosing  a  matrix  such  that 

(y(^).y)rFF'(y^^^y)  >  1,  ^»l,...,2n 

y^y  -  0  =— >  (Fy)'y  “  0  (39.2) 

S  (y^"^^ -y) 'FF' (y^"^^ -y) -1  - min 

i-1 

Note  that  0  may  happen. 

b)  Spherical  approximations:  An  outer  spherical  approximation 

j^(sph,2)^j^  is  defined,  cf.  (38.1a),  by 
o  o 

{yeE^;  y'y>0,  |  |y-y|  |  <  P  j  ) .  (^0) 

j*»2,  where  the  radius  given  by 

p  max  ||y^'*^-y||  (40.1) 

l<i<2n 

and  j|.||  denotes  the  Euclidean  norm.  Likewise,  an  "inner”  spherical  ap¬ 
proximation  of  can  be  determined,  cf.  (40),  by  choosing  the 

radius 

p  :»  min  1  -y  I  1  ,  (40.2) 

^  l<i<2n 

where,  as  for  the  relation  \  K  fpf  is  not  excluded  in 

o  o'  o'^^  .  h 

general.  Of  course,  an  inner  spherical  approximation  ^  is  ob¬ 

tained  by  selecting  the  radius 

min{  I  |y-y|  I  :  yeSK^),  (40.3) 

where  dK^  denotes  the  boundary  of  K^. 

Having  an  approximation  of  as  described  by  (39)  and  (40),  the 
cone  can  be  approximated  now  by  the  cone 

Y^:-  {Ay:  A>0,  yeK^). 


(41) 


tions  for<  y. 

Having  an  approximation  of  Y^,  the  probability  function  P-“P(A,X) 
given  by  (33.1)  can  be  approximated  then,  cf.  (36.2),  by 


P(A,X):-  P 


CF  (w,X)-AR^(tj) 
F^(w.X)-F^((o.X) 


Hence,  in  the  above  case  from  (42) -(42. lb)  we  obtain 

P(A,X)  -  P(| IT' (I lyl l^y(w)-y'y(w)  y) | |  <  y(w)'y) , 

2-2  1/2 

P(A,X)  -  P(lly(a))||  <  - y(o))'y), 


(43.1a) 

(43.1b) 
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resp . ,  where 


CF^(w,X) 


AR  (w) 
o 


Finally,  we  have  to  determine  the 


f'^(w.X)  -  F^((J,X) 


erating  vector  yGY^,  y+0.  Having  y 


(43.2) 

axis  g“(Ay:  A>0)  by  selecting  a  gen- 
(k) 

2  cr.y'^  ^  with  a,^0,  k-l,...,2n, 
k-1  ^  ^ 


condition  (38.1)  can  be  fulfilled  by  choosing  the  coefficients  a»(a^,.. 

. )'  such  that 
2n 


Y^a  >  0, 


where  the  2nx2n  matrix  Y^ 

o 


is  given  by 


(44) 


(44.1) 


where  is  the  nxn  identity  matrix.  Thus,  we  may  select  then  a  such  that 

(Jl)  *  - 

min  y  y  is  maximized,  which  yields  the  linear  program 
l<i<2n 

max  t  (45) 


a>0,  t>0. 


(45.1) 

(45.2) 


where  a^  *=  ^^ol'  %2  *  '  *  *  *  ^o2n^  ^  ^  given  2n-vector  having  positive  com¬ 
ponents,  and  the  constraint  ot<a^  is  imposed  because  the  direction  of  y  is 
needed  only. 

On  the  other  hand,  for  any  vector  y,  y+0,  we  consider  the  hyperplane  E 
defined  by  (37.1)  and  the  points  y^  \  i-1 . 2n,  on  E  defined  by  (38): 
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Since  the  axis  g={Ay:  A>0)  should  pass  through  the  center  of  the  cone 
Y^,  the  deviations  between  y  and  the  points  y  ,  i-l,...,2n,  should  be 
well-balanced.  Thus,  y  is  choosen  such  that  the  quantity 

max  I|y-y^'*^||  -  max  -  ’  Hyll^  > 

l^Z<2n  l<i<2n  (y'‘  y) 

cf.  (38)^,  is  minimized,  which  yields,  see  (38.1),  the  optimization  problem 

min  max  |  |y-y^  |  |  (47) 

l<i<2n 

s .  t . 

y"^  ^  y  >  0,  i-l,...,2n  (47.1) 

lly||“i;  (47.2) 

since  only  the  direction  of  y  is  relevant  for  our  purposes,  the  norm  con¬ 
straint  I |y| |«1  is  added.  According  to  (46),  problem  (47)  is  equivalent  to 
the  program 

max  t  (47)' 

s .  t . 

(i)'- 

t  -  y  :g  0,  i=l,...,2n  (47.1)' 

wy^’w 

t2t0,  yGY^.  ||y||-l. 


(47.2)' 
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Note  that 

(Jt)  - 

^  (pj  -  cos  \y), 

(i) 

is  the  cosine  of  the  angle  between  the  vectors  y  and  y. 

(k) 

Representing  y  e  by  y  -  S  y^  ^ k“l,...,2n,  we  observe 

k-l 

that  (47)'  is  closely  related  to  (45). 


8.  Sensitivity  analysis  of  probabilities  of  survival/failure 

In  the  following  we  consider  the  sensitivity  of  the  probability  of  sur¬ 
vival/failure  with  respect  to  the  variables  (A ,X, z) ,  i.e.,  with  respect  to 
the  nij^-vector  A  of  deterministic  load  factors,  the  r-vector  X  of  design 
variables,  and  the  (n-m) -vector  z  of  deterministic  redundants  in  the  n- 
vector  F“=F(ci;)  of  member  (element)  loads.  As  developed  in  [8-13],  deriva¬ 
tives  of  the  probability  functions  P“P(A,X),  P«=P(A,X,z)  can  be  obtained 
-  under  weak  mathematical  assumptions  -  by  means  of  the  Transformation 
Method  in  combination  with  a  stochastic  completion  technique. 

8.1.  The  probability  functions  (10.3)-(10.5) .  In  order  to  show  the  dif¬ 
ferentiation  of  the  probability  functions  given  by  (10 . 3) - (10 . 5)  it  is 
sufficient  to  consider  probability  functions  of  the  type 

P(A,X):=  P(AR^(w)  e  [A(X)  ,  A(X)a^(w)]),  (48) 

°  ^  L  U 

where  we  suppose  that  the  random  vectors  R^(w),  (a  (cj),a  (a?))  have  suffi¬ 
ciently  smooth  probability  densities  <p^(R^)  ,  ^V>(a^,a^)  on  re¬ 

spectively.  Hence,  using  also  (33) ,  we  have  that 


P(A,X)  »  P(A,A),  A=A(X)  »  (A^(X),...,A^(X))', 
cf.  (2.3),  where 

P(A,A):=  /  (p(R^)V'(ff^,a^)dR^d«^^dcr^ 

AR  eClA.a^.A.ff'^] 
o  ^  d  d  ^ 


(CA^a^-AR 
d  o 

A  U  .  L 

V 


/  (p(R^)^(cr^,ff'^)dR^dff^da^. 

G  Y 

o 


(48.1) 


(48.2) 


Applying,  for  given  variables  (A, A)  with  AfO,  A^>0,  i~l , . . . ,n«n^ , 
resp.,  to  (48.2)  the  transformation  T.  (R,F^,F^)  - (R  ,a^,a^)  in 

9ti  o 

X  X  R  defined  by 


R 

o 


L 

a  :  = 


u 


A'V 


(48.3) 


we  obtain 


(48.4) 


P(A,A)  -  / 


ReC[F^,F^] 
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,?(iR)^(A’^F^.A'^F^) 


L  U 

—  n  dRdF  dF  . 

.m  ^  ,  .2 


\xr  j“l  A, 


Under  weak  assumptions  [10], [13]  we  may  interchange  in  (48.4)  diffe¬ 
rentiation  and  integration,  hence. 


d? 

dX 


(X,A)-  -  X  /  ni<p(|))^(A'V,A'V) 

ReC[F^’F^] 

1  1  T  T1 

X  n  dRdF  dF  . 

m  ^  ,  2 


(49) 


|A|  j-1  A 

-1  •' 

Using  the  inverse  transformation  ,  (49)  yields 

||(A,A)-  -  X  •  /  divU^,p(R^))^(«7^,a^)dR^da^da^ 

AR  e  C[AjO-^,A^a^] 
o  ^  d  d 

-  -  i/  div(R^<p(R^))P(A,A|R^)dR^. 
with  the  conditional  probability  function 

P(A,A|R^):-  P(AR^  €  C[A^aV)  .A^a’^(w)  ] ) ,  R^  6 

ap 


(49.1) 


(49.2) 


Considering  now  the  derivative  j^(X,X)  with  respect  to  a  design  vari- 

K 


able  X^,  l<k<r,  by  means  of  (48.1)  we  find 
5P  ^  5P 


(50) 


j-1  -j  -k 

and,  using  again  (48.4),  by  interchanging  differentiation  and  integration 

there  we  get 

d?  ..  r 


dA. 

J 


A.  ^  ,  ^(f)t2*<A;V,A;lF'’) 

^  ReC[F^,F^] 


(51) 


,  ^(A^lp'-.A.lF'')^  .  2i!-(A.V,a;1f'’)^1 

So.  j  3<,^  J 

-  n  -  T  u 

X  n  ^  dRdF  dF  . 

1^1  J-1  A 

-1 

Using  the  inverse  T.  .  of  (48.3)  again,  from  (50.1)  we  get 

» A; 


dP 

dA 


;(^.A)  -  -  ^  f 

^  j  AR  EC[A,ff^,A,ff^] 


.p(R^)[2i^(<7^,a^) 


(51.1) 


,  ,  L  U,  L  .  a^  ,  L  U,  ,  L,  U 

+  -^(<r  ,o  )a  +  ,o  )cr  ]dR  d<7  Aa 


aa.‘ 

J 


j  aa.U' 


1  r  L  U.  (j)  ,,  L  ,  L  U.  .  L.  U 

”  ‘  A~  J  ^  ))P(A,A|ff  ,a  )da  da  , 

j 

with  the  conditional  probability  function 
P(A,A|a^,a^)  P(AR^(«)€C[A^o'^,A^a'^] ) 


(51.2) 


and 
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(0 . O.aJ'.O . O.a^.O . 0) ' ,  (5] 

L  U  ^ 

where  a.  , a.  are  placed  at  the  j-th,  (n-fj)-th  position,  respectively. 

J  J  L  U 

According  to  the  representation  (31)  of  the  event  (ReC[F  ,F  ]),  the 

conditional  probability  functions  P(A,A|R^),  P(A , A| can  be  repre¬ 
sented,  cf.  (31) -(33.1), (48.2),  by 
P(A,A|R^)  »  P(y(c.|R^)  e  Y  ), 

o  o  u  ^ 

P(A,A|cr'^,er^)  -  P(y(w| £t^,<7^)  e  Y^) , 
where 


(51.3) 


y(-|R„)H  ^  u 


A^cr  (w))  -  A^ff  (w) 


,  ,  L  U,  f  "V  - 
y^c  .a  ):  = 

\  V  -  V 


(52.1) 


Hence,  corresponding  to  the  approximation  P(A,A)  of  P(A,A),  cf.  (43)- 
(43.2),  the  above  conditional  probability  functions  (52)  can  be  approxi¬ 
mated  by 


P(A,A|R^):=  P(y(ulR^)  €  Y^), 
P(A,A|cr^,cr^)  :=  P(y(to|a^,a^)  G  Y^)  , 


(52.2) 


where  Y^  is  the  approximation  of  the  convex  cone  as  described  in  Section  7. 

By  a  slight  modification  in  the  equations  (49 . 1) , (51 . 1) ,  resp.,  the  de¬ 
rivatives  of  P=P(A,A)  can  also  be  represented  by  means  of  expectations: 
Theorem  8.1.  Under  the  assumptions  in  Section  8  we  get 
op  T  div(R  (a))(p(R  (w))) 

-  -  X  ^ ,  A  ^ 

o'  (w),A^ff  (w)]  ^^2) 

div(R  (w)(p(R  (w))) 

-  -  I  ^ - P(».MR<,<»)). 


1  ^  M  A'l  ^  T7  div((ff^(w)  ,ff^(w))^^  V(o^(w)  ,o^(w)))  ("iT  11 

b)  aXT^A.A)  -  -  ^  E  - L,  .  U,  ,, 

J  J  (w),(7  (W)) 

X  1  (AR  (a>)) 

C[A^a  (w),A^ff  (o))] 

.  1  E  diy.C(<T^((.).,aV))^^^^  ?(X,A\ah^),o^M) 


Remark  8.1. 

a)  Selecting  fixed  variables  (A, A)  such  that  AfO,  A.>0,  i«l,...,n,  and 

-1  ^  -1 
applying  then  -  instead  of  T.  ^  -  the  inverse  transformation  T  to 

(A, A) 


(49) ,  we  obtain 
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§1 

dX 


(A, A) 


; 


AR^e 


A  O  A  O  A  O 


d? 


//A-lx  L  a-1a  U,  ,A,m  !!  ,  L ,  U 

X  ^(A^  k^a  ,A^  k^a  )  |-(  II  (^)  dR^da  da  , 


3P 

dX’ 


and  ^(A,A)  can  be  represented  in  the  same  way.  Hence,  the  derivatives 
d? 

may  be  represented  by  integrals  over  the  fixed  domain  B 
j  (A,  A) 

{ (R^,cr^,<7^)  :  AR^  G  C[A^(7^,A^a^]  )  in  the  (R^ , a^, a^) -space . 

b)  Since  the  domain  of  integration  in  the  integral  representation  (49), 

5  P  5  P 

(51),  resp.,  of  — ,  is  independent  of  the  variables  A,  A»“(A^,  .  .  .  , A^) '  , 

the  higher  order  derivatives  -  of  arbitrary  order  -  of  P«P(A,A)  can  be  ob¬ 
tained  by  further  differentiation  of  the  equations  (49) , (51)  with  respect 
to  A,Aj ,l<j<n. 

c)  Having  the  mean  value  representations  (53), (53.1),  gradient  estimates 
-  as  well  as  estimates  of  the  probability  function  itself  -  can  be  obtained 
by  a  suitable  sampling  procedure. 

8.2  The  probability  function  (18.1),  Corresponding  to  (48)-(48.2)  we 
note  first  that 

P(X,z)  -  P(A,z),  A-A(X)  «  (A^(X),...,A^(X))^  (54) 

cf.  (2.3),  where 

r 


P(A,z):-  P 


C*^(R(w)-CjjZ)  S  Aj^^aj(w)l 


(54.1) 


Supposing  that  the  random  vectors  R(w) , (a^(w) ,a^(a)) )  have  sufficiently 


T  IT  m 

smooth  densities  <p==(p(R)  ,  ^«=^(a  ,a  )  on  31  ,31  ,  resp.  ,  and  defining  here 

the  integral  transformation 


>  (R,a^,a^)  by 


R:-  CjjZ  +  S 

L.  .  -1_L  L 

‘'l  ~  ^I.d  ^I’ 

U.  A  U 


(54.2) 


where  Aj>0,  l<j<n,  we  find 


I.d 


-1 


A  -V 

^I.d 


,^11.  d 


u 


(54.3) 


n 


X  n  ^dSAF^i-FKl- 

_  J 

Having  (54.3),  the  derivatives  of  P>-P(A,z)  -  of  various  orders  -  follow 
again  under  weak  assumptions  [10], [13]  by  interchanging  differentiation 
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and  integration  in  (54.3).  Hence,  corresponding  to  (51)"(51.3),  for 
j=l,2 . n  we  find 

J  j  A  /T^ 


T  T  IT  T  TT  ' 

X  P^(A,z|a  ,a  )d47  da  , 

where  P^-P^ (A, z | a^,a^)  is  the  conditional  probability  function  given  by 


P^(A,zla^,a^):-  <  Cj^(R(a>) -C^^z)  < 

Moreover,  for  jG{j^:  l<l<n-in)  we  get 

||-(A,z)  -  /V,?(R)'c  P(A,z|R)dR 

j  ^ 


(55.1) 


(55.2) 


1  r  ,a^  ,  L  U,  ^  ,  L  ,  L  U.  .  L.  U 

^  L  J  u  ^  ^  > 

i  A.,_  ,ff,  <z^.rx  jO'tt  3a.  aa. 

II, d  II  II, d  II  J  J 


(55.3) 


where 

Ik  ahw)  <  C:^(R-C  z)  <  A  ,a!J(«)\ 

P(A,z|R):-P  V  .  (55.3) 

\^II,d‘^II^“^  <  z  <  Ajj^^ajj(w)  / 

Remark  8.2. 

Using  the  inverse  transformation  T.-  ~)  ^  with  given  fixed  variables  z, 

QZ  ,  A 

A,  cf.  Remark  8.1,  we  also  obtain  integral  representations  of  the  deriva¬ 
tives  having  a  fixed  domain  of  integration  in  the  space  of  the  original 
(R,a^,a^) -variables . 


8.3.  The  probability  function  (26.1) 

According  to  the  definition  of  F  (a),X),F  (w,X)  given  in  Section  1  for 
different  cases,  the  probability  function  (26.1)  can  be  represented  by 


P(A,X)  -  P(a(w)'V(A,X)6  <  0  for  all  5eA^) , 


where 


Ro(a,) 


a(w)  :-l  a  (w)  I,  fi:-  u  I  , 
\a^(t<j)/  \  u  / 


(56.1) 


is  the  convex  polyhedron  of  elements  5  represented  by  the  system  of 
linear  equalities/inequalities  (23 . 1) - (23 . 3) ,  and 

V  =  V(A,X)  =  V(A,A(X),  W^CX),  W^(X),  wP(X))  (56.1) 
is  an  (m+2n)x(m+2n)  matrix  given  analytically;  in  case  of  trusses  we  have 
that 


V(A,X):- 


AI 

1 

-Ad(x) 

A„(X)  / 

(56.2) 


Having  the  extreme  points  6^  ,  i“l,...,i^,  of  A^,  we  also  get 

P(A,X)  -  P(a(w) 'V(A,X)5^'^^  <  0,  i=l . 


(56') 
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Since  the  number  2^  of  extreme  points  of  may  be  very  large,  and  the 
numerical  computation . of  is  very  time  consuming  in  gen¬ 

eral,  first  we  are  looking  for  upper  and  lower  bounds  of  P“P(A,X): 

Considering  an  arbitrary  sequence  of  elements 

^1*^2^  *  *  *  *  ^ j  *  *  *  *  * 

we  find  for  any  i/e  N  the  upper  bounds 

P^(A.X;5^,...,^^):»  P(a(a)) 'V(A,X)5j<0,  j=l,...,i/),  (57) 

where 

P^(A,X;5^,...,5^)  >  P^^^(A,X;5^,...,5^^^)  >  P(A,X).  (58) 

for  each  i/“l,2,,..  .  Obviously,  minimal  upper  bounds  P^(A,X)  are  obtained 
by  minimizing  (57)  with  respect  to  j-l,...,i/,  hence,  we  put 

P*(A,X):-  min{P^(A,X;5^,...,5^):  j-l,...,i/).  (59) 

By  means  of  optimum  upper  bounds  the  true  probability  P(A,X)  can  be 
reached  in  a  finite  number  of  steps: 

I^mma  8.1.  There  is  an  integer  ( A , X)  ,  such  that 

P*  (A,X)  P(A,X).  (60) 

o 

The  assertion  follows  from  the  inequalities 

P(A,X)  -  (A,x:5^^^ . 5  °) 

O 

<  P*(A,X)  <  P^(A,X;fi^^^...,«<''^ 

for  each  i/-l ,  2 ,  .  .  .  ,  i  . 

o 

According  to  [19],  suboptimal  upper  bounds  for  P(A,X)  can  be  obtained 
iteratively  as  follows: 

Stage  1.  Define 


P^(A,X):-  P^(A,X)  -  min{P^(A,X;5^):  (61a) 

"fc  it 

and  let  6--»5-(A,X)  denote  an  element  of  A  such  that 
11  o 

P*(A,X)  -  P^(A,X:5*(A,X)). 

Stage  I/.  Having  5^-5^  (A, X),  j-1 . t'-l,  for  v>\  define 

P*(A,X):-  min{P^(A,X:S*(A,X) . 5*_j^(A,X)  ;  SjEL^),  (61 

it  ^  if  St 

and  denote  by  6^^S^(\,X;S^ . an  optimal  solution  in  (61b).  Obvi¬ 

ously  we  have  that 

P*(A,X)  >  P*(A,X)  &  P(A,X),  »/-l,2 . 


(61b) 


Clearly,  a  big  advantage  in  (61b)  is  that  we  have  only  one  single  deci¬ 
sion  vector  6^,  whereas  in  (59)  we  have  to  deal  with  i/  decision  vectors 

On  the  other  hand,  with  the  suboptimal  upper  bounds  P*(A,X)  the 
exact  value  P(A,X)  can  not  be  reached  in  general  in  a  finite  number  of  steps. 
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According  to  (56')  we  find,  cf.  (29), (29.1), 

P(A,X)  -  1  -  P(F^U...UF^  ) 

o 

with  the  failure  domains  F^  given  by 

F^;-  (weO:  a(w) 'V(A,X)«^'^^>0) ,  Jl=l,  . . . 

Hence,  lower  bounds  for  P(A,X)  follow  by  applying  Bonferroni -bounds  [10] 

Jt 

o 

to  P(  ^  Fp).  These  bounds  can  be  described  by  means  of  an  estimate 
A  >e-i 

of  and  by  probability  functions  of  the  t3rpe 


Q^(A,X;6^,...,5^):-  P(a(a;) 'V(A  ,X)5  >0,  j-1,...,*/)  (62) 

similar  to  P^(A ,X; 5^ ,  .  .  . , 5^) ,  j-1,  . . .  ,i/. 

E.g.,  for  i/“l  we  get 
i 

®  (i) 

P(A,X)  >1-2  P(a(a))'V(A,X)5^  ^  >  0) 

i-1 
io 

-  1  -  S  Q  (A,X:5^'^^  >  1  -  i  Q*(A,X), 

-  -  i  O  J. 


where ,  for  , 2 , . . . , 

Q*(A,X):-  max{Q^(A,X;5^,...,5^):  l<j<i/)  . 

Consequently,  for  the  approximative  computation  of  the  probability 

P(A,X)  we  have  to  solve  optimization  problem  of  the  type 

minimize  P  (A ,X; 5- , . . . , ^  ). 

6.eA  ,l<j<j/  ^  ^ 

3  o* 


(62.1) 


(63) 


and 

maximize  Q  (A ,X; , . . . , 6  ).  (63.1) 

6  ,€A  ,l<j<i/  ^  ^ 

Moreover,  for  the  approximative  solution  of  reliability-oriented  opti¬ 
mization  problems  of  the  type 

max  P(A  ,X)  (64) 

XGD  ® 

with  a  given  load  factor  have  the  naxmln-problem 

maximize  minimize  P  (A  ,X;5i,...,5  ).  (64.1) 

XgD  5.eA  o  1 

J  o'  ^ 

Thus,  in  each  case  the  derivatives  of  the  probability  functions 
resp.,  are  needed. 

By  the  transformation 

V(A,X)6j,  j=l . »/,  (65) 

for  P^“P^(A,X;5^, . . . ,5^)  we  find  the  representation 

P^(A,X;5^ . S^)  -  W^(V(A,X)5j^ . V(A,X)fi^),  (66) 

where 

P(^ja(co)  <  0, 


(67) 
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and  can  be  represented  in  the  same  way.  Since  the  derivatives  of  V“V(A,X) 

can  be  obtained  analytically,  the  remaining  problem  is  the  differentiation 

of  W  . 
u 

8.3.1.  Exact  differentiation  formulas  in  case  of  1/  <  dim  a(.) 

In  many  practical  cases  the  stage  number  u  is  small,  hence,  the  assump¬ 
tion 

1/  <  d:-  dim  a(.)=  m+2n 

aw 

1/ 

is  not  too  restrictive.  For  the  computation  of  the  partial  derivative  -r-^ 

ou  . 


for  a  given  pair  (j,k),  l<j<i/,  l<k<d,  we  consider  a  partition  0 
of  the  i/Xd  matrix 
1 6  {\ 


jk 


0: 


(68) 


'2 

\KI 

such  that 

i)  0 ..  is  an  element  of  0_and 

jk  I 

ii)  rank  0^-1/. 

Supposing  then  that  the  random  vector  a*a(w)  has  a  probability  density 
f“f(a),  and  partitioning  the  vector  a€!R^  in  the  same  way  a^Ca^ja^^)  as  the 


(69) 


matrix  0,  by  the  integral  transformation 

.-1 


we  find 


W  ) 

j/  1  1/ 


; 

Pt+6ttPi 


SO 


fCe^^T-P 


r) 


<iPi 


I  |det0j.|  • 


(70) 


(71) 


I  Ill'll 

Obviously,  by  means  of  this  transformation,  the  domain  of  integration 
in  (71)  is  now  independent  of  the  element  and  -  of  course  -  also  in¬ 
dependent  of  all  other  elements  $  contained  in  0_ .  Thus,  if  the  density 

t  K  JL 

f  of  a(ct))  is  sufficiently  smooth  the  derivative  follows,  cf.  [10],  [13],  by 
interchanging  differentiation  and  integration: 

Theorem  8.2.  Under  appropriate  assumptions  [10], [13]  on  the  density 

Q 

f=f (a)  of  a(w) ,  the  partial  derivative  ,  is  given  by 


aw 


dB. 


<9 


.)• 


(a)  ■  (0'^'  )^aj^+f  (a)  (0’^’  )ji^)da. 


-W 


.-1' 


jk  0a<O 

which  can  be  represented  also  by  the  expectation 


aw 


30 


<0 


r 


Va  f(a(w))  ,  . ,  ,  - , 

■V-  -"«  f(aw)  > 


I  'jk'  [0a<O] 


-1'  -1' 
where  (0,  ).  denotes  the  i-th  row  of  0, 

I  J  ■’  I 


(72) 


(a(a))), 

(72.1) 


-  26  - 

aw 

Corollary  8.1.  The  derivatives  —  of  with  respect  to  elements  0 

LK 

of  have  the  same  form. 

8.3.2.  Approximative  derivatives  of 

If  condition  (71)  can  not  be  fulfilled,  e.g.  in  the  case  i/>d,  then  ap¬ 
proximative  derivatives  of  -  of  an  arbitrary  high  accuracy  -  can  be  ob¬ 
tained  by  the  following  stochastic  completion  technique  [8]  ,  [13]  : 

Let  z^-Zj(w),  denote  real  random  variables  such  that 

i)  Zj  (w) ,  j*=l,...,i/,  and  a(w)  are  stochastically  independent 

ii)  Ez.(w)«0,  and  z.(w)  has  a  continuous 

^  ^  (73) 

probability  density  \^j-0j(t),  j-*l,...,i/. 

By  means  of  the  "stochastic  completion  terms"  z^-Zj  (u>)  ,  j=l,.,.,j/,  for 

we  obtain  the  approximative  probability  function 


. . .  ,e^):=  P(«j'a(w)  +  ZjCw)  ^  0,  j-1 . v) . 


We  find 


. ^  -  ^j'a(w),  j-1, 


-E  n  ^A-e>a(o,)), 

j-1  J 

where  denotes  the  distribution  function  of  z . (w) :  if  z.(w)  has  the 
J  J  1 

normal  distribution  N(0,a.),  then 

^  1 

. E  n  (75.1) 

^  j-1  ‘^j  J 

where  ^  is  the  distribution  function  of  N(0,1). 

The  partial  derivatives  of  W^  read 


d  e,. 


n  t^(-^^'a(w))V'j  (-Sj 'a(w))aj^(w)  , 


j=l,...,j/,  k=l,...,d,  and  higher  order  derivatives  of  W^  can  be  obtained 
in  the  same  way.  If 

z(a))  =  (Zj^(w)  ,  .  .  .  ,2^(w))  '  - 0  w.p.l  (77 

(with  probability  one),  then  under  some  regularity  assumptions  [9], [13] 

. 


. V 


a  . 

Jk 
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